clear all
%clc

warning('off','MATLAB:dispatcher:UnresolvedFunctionHandle')

nMC = [1:250];

theta_all = [];
ratio_trimmed_mean_all_second = [];
ratio_trimmed_mean_all_first = [];
trimmed_mean_of_ratios_all = [];


for i=nMC
    try
    tempp = num2str(i);
    temp = strcat('MC_',tempp);
    temp2 = strcat(temp,'.mat');
    
    load(temp2)
    
    theta_all = [theta_all;
        theta_NPD'];
    
    ratio_trimmed_mean_all_second = [ratio_trimmed_mean_all_second;
                              ratio_trimmed_mean_second];

    ratio_trimmed_mean_all_first = [ratio_trimmed_mean_all_first;
                              ratio_trimmed_mean_first];

    trimmed_mean_of_ratios_all = [trimmed_mean_of_ratios_all;
                                  trimmed_mean_of_ratios];

    catch
    end

end

theta_NPD = theta_all';

display('*** RATIO OF TRIMMED MEANS -- SECOND ***')
mean(ratio_trimmed_mean_all_second)
std(ratio_trimmed_mean_all_second)

display('*** RATIO OF TRIMMED MEANS -- FIRST ***')
mean(ratio_trimmed_mean_all_first)
std(ratio_trimmed_mean_all_first)
